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SUMMARY 


I  Describe^hereir^are  the  method  of  obtaining  photographic  meteor  data  in  the  mass 
range  of  10"^o  10®  gram  and  the  treatment  of  these  data  to  provide  useful  engineering 
relations  for  protecting  large  space  vehicle  components  against  meteoroid  impach  The 
assumptions  needed  to  convert  these  data  to  meteoroid  masses  and  densities  are  dis¬ 
cussed,  and  the  effects  of  these  assumptions  on  design  r^fations  are  shown.  The  inferred 
nature  and  values  of  meteoroid  density,  mass,  and  velocity  are  described,  and  an  inter¬ 
nally  consistent  meteoroid  mass  flux  model  is  introduced  for  the  photographic  meteoroid 
range Earth  shielding  factors  as  a  function  of  satellite  altitude  are  presented  for  use 
with  the  proposed  model.  Values  of  average  meteoroid  density  and  average  meteoroid 
velocity  relative  to  the  Earth  at  the  Earth’s  atmosphere  lower  than  values  previously 
used  are  indicated  by  an  analysis  of  the  available  photographic  meteor  data.  Photo¬ 
graphic  meteor  astronomy  is  shown  to  be  a  significant  source  of  data  in  describing  the 
meteoroid  hazard  in  space.!  . 
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INTRODUCTION 


The  space  vehicles  and  power  generation  systems  of  the  future  can  be  expected  to 
have  stringent  requirements  of  long  duration,  large  surface  areas,  and  high  probability 
of  success  placed  upon  them.  Components  of  these  flight  systems,  especially  those  con¬ 
taining  fluid  circuits  and  life  support  areas,  will  have  to  be  protected  against  critical 
damage  due  to  meteoroid  impact.  Components  with  especially  large  exposed  surface 
areas  are  the  propellant  tanks  of  large  chemical  and  nuclear  rockets  and  the  waste-heat 
radiators  of  power  generation  and  coolant  systems. 

The  critical  meteoroid  mass,  or  the  mass  of  the  largest  meteoroid  that  the  vehicle 
component  must  be  protected  against,  increases  with  the  vulnerable  area,  the  exposure 


time,  and  the  probability  of  survival  without  critical  damage  (e.  g. ,  ref.  1).  This  rela¬ 
tion  indicates  that  large  space  system  components  must  be  protected  against  meteoroid 
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particles  with  masses  of  the  order  of  10  to  10  gram.  Of  particular  interest  to  the 
designer  of  such  systems  are  meteoroid  influx  rate,  density,  structure,  impact  velocity, 
and  direction  of  flight  in  this  particular  mass  range.  Since  the  design  engineer  cannot 
predict  the  exact  conditions  of  these  variables,  he  must  rely  on  the  probability  of  meeting 
specified  conditions. 

Numerous  satellites  and  rockets  have  been  launched  in  an  attempt  to  count  directly 

the  particle  flux  in  the  near-Earth  environment  (e.  g. ,  ref.  2).  The  recent  experiment. 

Explorer  XVI,  was  the  first  satellite  to  measure  directly  perforation  rates  in  pressurized 

cells.  The  energy- sensitive  detectors  of  Explorer  XVI  recorded  particle  flux  rates  two 

to  three  orders  of  magnitude  lower  than  the  rates  obtained  by  previously  flown  momentum- 
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sensitive  experiments  for  the  10  to  10  gram  range  (ref.  3).  Preliminary  reduction  of 

these  data  in  terms  of  a  geocentric  and  a  heliocentric  distribution  of  near-Earth  particles 

indicates  a  possible  reduction  in  the  apparent  discrepancy. 
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Unfortunately,  satellite  measurements  in  the  10  to  10  gram  range  are  from  four 
to  seven  orders  of  magnitude  from  the  range  of  interest  for  large  space  vehicle  compo¬ 
nents,  and  the  data  cannot  justifiably  be  extrapolated  into  this  region.  An  experiment  of 
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the  Explorer  XVI  type  would  be  highly  desirable  in  the  10  to  10  gram  range.  Accord¬ 
ing  to  the  best  currently  available  estimates  of  the  flux,  however,  a  satellite  with  approx- 

5 

imately  10  square  feet  of  vulnerable  area  and  an  exposure  time  of  approximately  1  year 
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would  yield  only  one  impact  of  a  10  gram  particle.  Directly  applicable  satellite  experi¬ 
ments  thus  appear  to  be  highly  improbable  at  the  moment. 

The  design  engineer  must  therefore  rely  on  indirect  meteoroid  data  obtained  by  as¬ 
tronomical  means  (photographic,  visual,  and  radar),  which  measure  particles  with 

_4 

masses  10  gram  and  larger.  These  meteor  data  fall  well  within  the  range  of  interest 
for  the  large  long-duration  space  vehicles,  and  no  unjustifiable  extrapolation  is  neces¬ 
sary.  Meteor  astronomy  is  the  general  term  defining  the  study  of  meteors.  Specifically, 
a  meteor  is  the  streak  of  light  caused  by  a  meteoroid  traveling  at  high  velocity  through 
the  Earth's  atmosphere.  The  study  of  meteor  astronomy  is  therefore  specifically  limited 
(by  definition)  to  these  large  particles  that  leave  visible  evidence  of  their  passage.  Gen¬ 
eral  usage  of  the  term,  however,  has  led  to  the  incorporation  of  all  particles  (micro¬ 
meteoroids  to  bolides)  under  this  general  heading;  however,  the  discussion  herein  will  be 
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limited  to  photographic -size  meteoroids  (10  g  and  larger). 

Considerable  meteor  data  are  available  in  the  literature,  and  various  estimates  of 
influx  rates,  densities,  and  velocities  based  on  these  data  have  been  proposed  (e.  g. , 
refs.  4  to  6).  Radical  differences  in  influx  rates  and  densities  of  meteoroids  from  pro¬ 
posal  to  proposal  lead  one  to  be  concerned  with  the  significance  and  accuracy  of  astro¬ 
nomically  derived  predictions.  Order  of  magnitude  differences  appear  in  these  various 
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proposals,  and  the  design  engineer  is  at  a  loss  to  understand  the  reason  for  such  differ¬ 
ences  and  the  shortcoming  of  proposing  average  density,  velocity,  and  influx  rate  as  uni¬ 
versally  applicable.  Use  of  these  various  proposals  for  space  vehicle  design  criteria  re¬ 
sults  in  large  differences  in  protection  requirements  and  weight  among  design  studies.  It 
is  important,  then,  to  understand  the  capabilities  and  limitations  of  astronomical  reduc¬ 
tion  procedures  and  how  to  interpret  the  data  presented  in  terms  of  the  particle  charac¬ 
teristics  and  flux. 

One  of  the  objectives  of  this  report  is  to  explain  how  photographic  meteor  data  have 
been  transformed  into  meteoroid  densities  and  influx  rates  through  the  laws  of  meteor 
physics.  The  assumptions  necessary  in  the  reduction  of  the  data  are  discussed,  and  the 
confidence  that  can  be  placed  on  calculated  protection  requirements  is  shown.  In  addi¬ 
tion,  application  of  the  reduced  data  is  discussed  and  a  consistent  meteoroid  model  in  the 
photographic  region  is  presented. 


METEOR  PHYSICS 
Description 

Meteoritic  particles  are  collectively  referred  to  as  meteoroids.  Present  knowledge 
of  meteoroids  suggests  that  they  be  classified  into  two  major  groups,  with  their  origins 
designated  as  either  cometary  or  asteroidal.  The  majority  of  meteoroids  (in  excess  of 
90  percent)  are  considered  to  be  of  cometary  origin.  These  particles  represent  the  de¬ 
bris  left  in  space  after  cometary  passage.  The  mechanism  of  injection  and  the  type  of 
particle  injected  are  discussed  in  a  following  section.  At  this  point,  however,  it  is  well 
to  define  terms  associated  with  cometary  meteoroids.  A  meteoroid  shower  is  a  periodic 
increase  of  meteoroid  flux  and  represents  a  group  of  particles  traveling  in  a  definite, 
predictable  orbit  around  the  Sun.  The  fact  that  many  major  showers  have  been  corre¬ 
lated  with  known  comets  confirms  their  cometary  origin.  Sporadic  meteors  are  not  in 
predictable  orbits  and  represent  cometary  meteors  that  have  had  time  to  depart  appre¬ 
ciably  from  the  initial  cometary  orbits  because  of  perturbing  forces. 

Asteroidal  meteoroids  are  considered  to  be  of  different  origin  from  the  cometary 
meteoroids.  Their  origin  is  believed  to  be  from  the  asteroids,  a  number  of  large 
masses  (up  to  a  800-km  diam)  in  orbit  between  Mars  and  Jupiter.  As  shown  in  a  follow¬ 
ing  section,  the  asteroidal  meteoroids  are  of  a  much  different  structure  from  the  typical 
cometary  meteoroids.  Their  density  is  about  an  order  of  magnitude  greater  than  that  of 
the  cometary  meteoroids,  and  sometimes  they  are  found  to  be  90  percent  iron.  In  gen¬ 
eral,  they  are  a  compact  stony  type  of  particle,  and  frequently  they  survive  the  passage 
through  the  Earth's  atmosphere  and  are  found  on  the  Earth's  surface.  Data  are  obtained 
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on  the  larger  cometary  (and  very  infrequently  on  the  asteroidal)  meteors  by  photographing 
the  meteor  trail. 

Ideally,  photographic  meteor  data  supply  four  types  of  information:  particle  velocity, 
deceleration,  trajectory,  and  photographic  magnitude.  These  data  are  obtained  by  photo¬ 
graphs  of  the  meteor  taken  simultaneously  with  two  Baker  Super- Schmidt  cameras  for 
which  the  base  line  distance  between  them  is  accurately  known  (approx.  40  km).  In  each 
camera  is  a  rotating  shutter  with  a  known  period,  and  the  area  of  the  sky  photographed  by 
each  camera  is  known  (typically,  5980  sq  km  at  a  90-km  altitude) .  Unfortunately,  not  all 
the  preceding  information  can  be  obtained  for  all  meteors  photographed  with  this  equip¬ 
ment.  The  fainter  meteors  and  those  that  do  not  have  a  complete  trail  on  at  least  one  of 
the  plates  cause  complications  in  the  determination  of  the  particle  characteristics.  The 
fainter  meteors  have  very  short  visible  trails,  which  makes  particle  deceleration  mea¬ 
surements  virtually  impossible. 

From  the  photographic  plates  of  a  meteor  that  is  bright  enough  and  the  trail  of  which 
is  long  enough,  the  velocity  can  be  calculated  by  measuring  the  distance  between  breaks 
in  the  meteor  trail  caused  by  the  rotating  shutter.  The  deceleration  is  found  by  measur¬ 
ing  the  change  in  these  distances  between  breaks,  and  the  photographic  magnitude  is 
found  by  comparing  the  meteor  brightness  with  the  star  background.  Finally,  each  cam¬ 
era  and  the  line  of  the  meteor  determine  a  plane  in  space  containing  the  meteor  path. 

The  intersection  of  the  two  planes  represents  the  meteor  path  in  space.  Extension  of  the 
planes  allows  the  determination  of  the  meteor  trajectory  (i.  e. ,  height  and  radiant) 

(ref.  7). 


Basic  Relations 

An  understanding  of  the  physical  equations  governing  the  theoretical  meteoric  pro¬ 
cess  is  desirable  in  order  to  understand  the  procedure  and  assumptions  made  in  calcu¬ 
lating  meteor  densities  and  influx  rates  from  the  basic  meteor  data.  Basically,  these 
equations  include  conservation  of  momentum,  conservation  of  energy,  luminosity  rela¬ 
tion,  and  various  combinations  of  these  equations  (ref.  8).  In  order  to  calculate  meteor¬ 
oid  densities  and  influx  rates,  however,  use  of  only  the  conservation  of  momentum  and 
the  luminosity  relation  are  required.  The  equations  presented  are  valid  only  for  solid 
compact  bodies  traveling  through  the  atmosphere.  Fragmentation,  or  the  breakup  of  the 
meteoroid  as  it  passes  through  the  atmosphere,  complicates  the  use  of  these  equations 
as  shown  in  the  following  section. 

Conservation  of  momentum.  -  Suppose  that  a  compact  meteoroid  of  mass  m  and 
density  enters  the  atmosphere  with  velocity  v.  (All  symbols  are  defined  in  the  ap¬ 
pendix.  )  Air  particles  are  trapped  on  or  near  the  surface  of  the  meteoroid,  impart  to  it 
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the  energy  of  collision,  and  decrease  its  forward  momentum.  The  air  particles  are 
given  a  forward  velocity  while  the  meteoroid  suffers  a  deceleration  dv/dt;  thus,  by  the 
principle  of  conservation  of  momentum, 

m  —  =  -  riSiov^  (1) 

dt 


where  p  is  the  density  of  the  atmosphere,  r  is  the  particle  drag  coefficient,  and  S  is 
the  effective  cross-sectional  area  of  the  meteoroid.  One  of  the  variables  in  equation  (1) 
can  be  eliminated  by  expressing  the  effective  cross-sectional  area  S  as  a  function  of  the 
mass  m.  An  assumption,  however,  must  be  made  as  to  the  shape  of  the  meteoroid.  If, 
for  example,  the  meteoroid  is  a  sphere  of  density  p^. 


(2) 

1/3 

In  general,  the  meteoroid  is  an  irregularly  shaped  object,  and  the  term  (97r/16)  '  in 
equation  (2)  is  replaced  by  a  dimensionless  shape  factor  A,  which  differs  from  body  to 
body.  With  these  transformations,  equation  (1)  can  be  written  as  the  first  fundamental 
equation  in  meteor  theory,  called  the  drag  equation; 


dv 

dt 


=  -rAp;;?/2m-V3p,2 


(3) 


Luminosity  relation.  -  The  blackbody  radiation  from  the  surface  of  the  meteoroid 
represents  only  a  small  fraction  of  the  luminosity  for  the  relatively  small  masses  of  the 
photographic  or  visual  meteors  under  consideration.  The  principal  source  of  the  visual 
or  photographic  radiation  occurs  in  a  coma  of  vaporized  meteoric  material  around  and 
behind  the  meteoroid  itself.  Practically  all  the  original  energy  is  expended  as  the  cloud 
of  meteoric  atoms  collides  with  molecules  of  the  atmosphere  after  the  meteoric  material 
has  left  the  main  body.  The  luminosity  I  is  therefore  assumed  to  be  equal  to  the  energy 
of  the  mass  lost  per  second  multiplied  by  a  luminous  efficiency  factor  t.  The  luminosity 
is  then  expressed 


I 


1  ^v2t 

2  dt 


(4) 


Thus,  there  are  two  basic  equations,  (3)  and  (4),  with  mass,  velocity,  and  time  as  three 
variables;  there  are  four  unknown  quantities:  P,  p^,  t,  and  A.  Two  of  the  variables. 
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V  and  t,  can  be  determined  from  photographic  data,  and  the  mass  m  is  still  to  be  found. 
The  integration  of  equation  (4)  yields 


The  terminal  mass  of  the  meteor  can  be  shown  to  be  nearly  zero,  so  that 


m  = 


(6) 


As  shown  in  a  following  section,  the  luminous  efficiency  t  is  related  to  a  constant  t^, 
called  the  luminosity  coefficient,  by  the  relation  r  =  t^v.  Then, 


m  = 


2 


(7) 


The  mass  referred  to  in  equation  (7)  is  the  instantaneous  mass  at  time  t  (the  lower  limit 
of  the  integral) . 

Substituting  the  mass  from  equation  (7)  into  equation  (3)  and  rearranging  terms 
yields 


T 


O 


P 


2 

m 


(8) 


Since  the  mass  in  equation  (3)  is  an  instantaneous  value,  this  equation  is  valid  at  a  point 
corresponding  to  the  lower  limit  of  the  integral.  The  values  of  v,  dv/dt,  and  p  must 
be  the  values  at  this  point  corresponding  to  time  t. 

With  these  limitations,  the  right  side  of  equation  (8)  can  be  determined  from  photo¬ 
graphic  data.  The  velocity,  light  intensity,  deceleration,  duration  of  the  meteor,  and 
atmospheric  density  (determined  from  the  known  particle  height)  are  all  measured  quan¬ 
tities.  In  order  to  determine  the  density  of  the  meteoroid,  t  and  TA  must  thus  be  de¬ 
termined  (or,  conversely,  if  and  TA  are  known,  the  luminous  efficiency  for  a  par¬ 
ticular  velocity  can  be  determined).  The  preceding  discussion  indicates  that  meteoroid 
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densities  can  be  calculated  readily  from  equation  (8)  from  photographic  data  along  with 
estimated  values  of  the  luminous  efficiency  and  the  drag- area  coefficient,  as  long  as  the 
meteoroids  remain  compact  single  bodies. 

Fragmentation.  -  The  typical  meteoroid  does  not  meet  the  conditions  of  compact  body 
flight  through  the  atmosphere.  The  development  of  a  "correction  factor"  to  be  applied  to 
the  equations  of  compact  meteor  flight  was  therefore  necessary.  This  correction  factor, 
called  the  fragmentation  index,  was  developed  by  Jacchia  (ref.  9)  as  an  answer  to  the 
"faint  meteor  anomaly.  "  This  anomaly  was  discovered  when  data  from  the  more  recent 
Baker  Super-Schmidt  cameras  were  compared  with  the  predictions  of  the  "ballistic" 
theory  just  presented.  (The  Baker  Super- Schmidt  cameras  are  capable  of  photographing 
much  smaller,  i.  e. ,  fainter,  meteors  than  cameras  previously  used. )  The  durations  of 
these  meteors  were  shorter,  and  their  decelerations  increased  more  rapidly  than  pre¬ 
dicted  by  theory. 

These  and  other  observed  anomalies  can  be  explained  by  the  process  of  progressive 
fragmentation  proposed  by  Jacchia.  It  is  theorized  that  the  meteor  fragments  progres¬ 
sively  along  its  trail.  Small  fragments  can  be  detached  from  the  surface  of  larger  me¬ 
teor  bodies  without  disrupting  their  unity,  but  if  fragments  of  similar  size  are  detached 
from  small  bodies,  complete  disruption  into  a  cluster  of  fragments  results.  This  phe¬ 
nomenon  explains  the  higher  frequency  of  occurrence  of  anomalies  in  the  fainter  meteor 
trails. 

Jacchia  developed  a  quantity,  the  fragmentation  index,  as  a  convenient  measure  of 
the  progressive  departure  of  the  observed  meteor  deceleration  from  single-body  theory. 
The  fragmentation  index  y  is  empirically  determined  and  is  defined  as 

^2^  (9) 

^t  / 

where  a^^g  is  the  observed  acceleration,  a^  is  the  acceleration  computed  from  the 
drag  equation  (3),  and  s  is  the  mass  loss  parameter,  defined  as 

s  =  log 

where  m^  is  the  original  mass  of  the  meteor  before  it  entered  the  Earth's  atmosphere, 
and  m  is  the  instantaneous  mass  of  the  meteor.  Unfortunately,  y  is  not  easy  to  deter¬ 
mine,  because  it  involves  the  second  time  derivative  of  the  velocity;  however,  x  is  ap¬ 
proximately  constant  for  each  meteor  during  the  detectable  part  of  its  flight  (ref.  10). 

The  effect  of  fragmentation  on  equation  (8)  can  be  readily  shown  by  rewriting  equa¬ 
tion  (8)  as 


This  equation  was  derived,  as  explained  previously,  by  substituting  a  value  for  m  into 
the  drag  equation.  The  drag  equation  gives  the  theoretical  value  for  dv/dt;  equations  (8) 
and  (11)  are  thus  ideal  or  theoretical  equations.  In  this  idealized  case,  remains 

a  constant  along  the  meteor’s  path  (ref.  10),  and,  therefore,  as  long  as  consistent  values 
of  V,  p,  and  dv/dt  are  used,  the  calculated  values  of  should  remain  the  same. 

This  implies  that  the  observed  value  of  acceleration  is  equal  to  the  theoretical 

value  of  acceleration  a^  (as  given  by  the  drag  equation),  and  that  the  fragmentation  index 
X  as  seen  from  equation  (9)  is  zero.  In  reality,  equation  (11)  is 


2 

If  the  value  of  TApv  from  equation  (3)  is  substituted,  equation  (12)  becomes 


If  m  is  substituted  in  the  idealized  situation  from  equation  (7), 


Equation  (14)  further  shows  that  if  the  observed  acceleration  equals  the  theoretical  accel¬ 
eration,  (jq/p^^  equals  which  is  a  constant. 

Differentiating  the  logarithms  of  equation  (14)  yields 
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is  a  constant  for  a  given  meteor.  Substituting  equation  (9)  into  (15) 


since 

yields 


ideal 


(16) 


and  integrating  from  s^  to  s  yields 


Equation  (17)  shows  how  two  values  of  log  meteor  trajectory  correspond¬ 

ing  to  two  values  of  the  remaining  meteoroid  mass  differ  from  one  another.  Again,  in  the 
absence  of  fragmentation  (y  =  0),  the  two  values  are  the  same. 

The  calculated  value  of  log  ^ny  meteor  can  now  be  related  to  some  ref¬ 

erence  point  corresponding  to  s^.  The  data,  to  be  most  significant,  should  be  reduced 
to  the  value  of  s^  corresponding  to  the  beginning  of  fragmentation.  Unfortunately,  this 
value  is  not  generally  known,  and  the  data  are  corrected  to  a  value  of  s  near  the  begin¬ 
ning  of  the  meteor.  It  is  indicated  intuitively,  that  the  fragmentation  index  is  also  an  in¬ 
dex  of  credibility;  that  is,  the  data  from  meteors  with  very  little  fragmentation  are 
"corrected”  by  the  least  amount  and  represent  data  that  deviate  the  least  from  theoreti¬ 
cal  considerations. 

An  excellent  example  of  the  use  of  photographic  meteor  data  with  the  fragmentation 
correction  is  given  by  Verniani  (ref.  10).  Also  discussed  is  a  system  of  "weights"  by 
which  the  calculated  values  of  log  different  values  of  s  along  the  meteor 

trajectory  are  deemed  credible  or  not. 

With  the  introduction  of  the  fragmentation  correction,  it  is  now  generally  possible  to 
use  equation  (11)  to  determine  meteoroid  density  if  the  luminous  efficiency  t  and  the 
drag-area  coefficient  TA  are  known. 

Luminous  efficiency.  -  The  numerical  value  of  the  luminous  efficiency  t  is  perhaps 
the  largest  single  influencing  factor  in  the  reduction  of  meteor  data.  Estimates  of  t  ap¬ 
pearing  in  the  literature  have  differed  from  one  another  by  as  much  as  two  orders  of  mag¬ 
nitude.  Opik  (ref.  11)  has  offered  a  theoretical  evaluation  of  t  for  both  solid  and  dust- 
ball  photographic  meteors.  His  analysis  is  lengthy  and  involved,  and  by  far  the  most 
complete  work  to  date.  The  analysis  calculates  the  radiation  energy  emitted  in  the  wave¬ 
length  range  4500  to  5700  A  by  iron  atoms  in  a  cloud  through  which  a  stream  of  nitrogen 
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Figure  1.  -  Variation  of  theoretical  luminous  efficiency  with  velocity  according  to  Opik 
(ref.  11). 


atoms  is  moving  with  a  series  of  specific  high  velocities.  For  the  diluted  coma  condition 
(dust  ball),  he  assumed  that  the  physical  process  is  better  represented  by  single  iron 
atoms  moving  through  a  cloud  of  nitrogen  atoms,  which  he  shows  produces  less  total  ra¬ 
diation  energy. 

The  results  of  Opik's  work  are  shown  in  figure  1,  in  which  luminous  efficiency  is 
plotted  as  a  function  of  meteor  velocity  relative  to  the  Earth.  Note  the  radical  depen¬ 
dence  of  T  on  velocity  at  values  less  than  15  kilometers  per  second.  Opik’s  work  pre¬ 
dicts  a  value  of  zero  for  luminous  efficiency  at  velocities  less  than  5.  2  kilometers  per 
second.  At  high  velocities,  however,  the  luminous  efficiency  is  very  nearly  constant  for 
both  the  solid  particle  and  the  dust  ball.  Opik's  theoretical  values  have  been  used  quite 
extensively  by  investigators  in  the  field  of  meteor  astronomy  (e.  g. ,  refs.  12  and  13). 

An  attempt  has  recently  been  made  to  measure  the  luminous  efficiency  of  an  artifi¬ 
cial  iron  meteor  of  known  mass  and  density  (ref.  14).  In  this  experiment,  a  seven-stage 
rocket  (Trailblazer  I)  was  employed  to  accelerate  a  stainless- steel  projectile  of  approxi¬ 
mately  2.  2  grams  into  the  atmosphere  at  10  kilometers  per  second.  The  entry  of  the 
projectile  was  recorded  photographically,  and  methods  of  photographic  reduction,  well 
documented  (e.  g. ,  ref.  7),  were  employed  to  obtain  the  luminous  efficiency  for  iron  me¬ 
teors.  This  value  was  further  reduced  to  the  efficiency  of  stony  particles  by  assuming 
that  the  stony  particle  was  15  percent  iron  and  that  all  the  visible  light  given  off  was  due 
to  iron  (ref.  11).  This  value  is  then  assumed  typical  of  meteoroids  that  are,  for  the  lack 
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of  better  information,  assumed  stony  in  composition. 

As  another  source  of  data,  values  of  luminous  efficiency  for  three  photographed  as- 
teroidal  meteors  were  calculated  by  Cook,  Jacchia,  and  McCrosky  (ref.  15).  Two  of  the 
meteors  were  assumed  to  be  stone.  The  other  was  identified  to  be  pure  iron  from  its 
available  spectrum.  Because  of  the  assumed  compact  nature  of  the  asteroidal  particles 
(as  opposed  to  the  » ’fluffy”  cometary  particles),  a  good  approximation  to  their  density 
was  made  and  the  luminous  efficiency  was  calculated  with  a  fair  degree  of  confidence. 

The  results  of  these  calculations  for  the  asteroidal  particles  are  shown  in  figure  2  along 
with  the  value  obtained  from  the  Trailblazer  experiment.  Fairly  good  agreement  with 
Opik’s  theoretical  calculation  for  a  solid  particle  is  seen.  Unfortunately,  however,  all 
these  experimental  points  are  in  the  low- velocity  range  of  the  curve  and  do  not  prove,  or 
disprove,  the  constancy  of  t  at  high  velocities  (>20  km/sec)  for  this  type  of  particle. 

Most  investigators  have  recognized  the  dependence  of  the  luminous  efficiency  of  a 
cometary  meteor  upon  its  velocity.  In  general,  for  analytical  purposes,  an  equation  of 
the  form 

T  =  TqV^  (1^) 

is  used  to  describe  this  dependence,  where  t  is  the  luminous  efficiency,  is  a  con¬ 
stant  (luminosity  coefficient),  and  v  is  the  meteoroid  velocity.  The  exponent  n  is  be¬ 
lieved  to  be  somewhat  less  than  one  (ref.  16).  For  simplicity  in  calculation,  however, 
and  because  of  the  fact  that  the  experimental  data  can  reasonably  be  fitted  within  its 
range  of  error,  most  investigators  use  1  for  the  exponent  n;  that  is, 

T  =  TqV  (19) 

Shown  also  on  figure  2  is  the  analytical  curve  for  t  according  to  equation  (19)  where 

T  =  10“^®-  (units  of  zero  magnitude  (visual)  g"^cm"^sec^),  given  by  Whipple  as  the 
o 

best  estimate  for  use  in  his  calculations  (ref.  17). 

Recently,  Verniani  (ref.  10)  made  an  extensive  evaluation  of  the  value  of  n  and  of 
the  absolute  value  of  for  cometary  meteors.  His  analysis  used  400  precisely  re¬ 
duced  meteors  and  represents  the  first  attempt  to  analyze  existing  data  on  cometary  me¬ 
teors  for  this  purpose.  His  analysis  takes  fragmentation  into  account  and  is  based  on  a 
knowledge  of  average  meteoroid  density.  Also,  as  indicated  in  the  next  section,  his 
analysis  treated  so  large  a  sampling  of  meteors  that  the  effects  of  unknown  individual  val¬ 
ues  of  the  drag- shape  factor  were  reduced.  Furthermore,  he  indicated  that  n  =  1  is  a 
good  representation  of  the  luminosity  law  over  the  complete  range  of  meteor  velocities 
(10  to  72  km/sec).  It  is  indicated  from  the  analysis  that  a  linear  relation  between  lumi¬ 
nous  efficiency  and  velocity  (eq.  (19))  may  be  a  better  representation  of  the  luminosity 
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law  for  cometary  meteors  than  that  proposed  by  Opik’s  dust-ball  curve.  The  analytical 
variation  for  t  obtained  from  Verniani  is  also  shown  in  figure  2. 

Drag- shape  factor  FA.  -  The  particle  drag  coefficient  r  has  been  shown  to  be  ap¬ 
proximately  equal  to  1  in  free  molecular  flow  and  approximately  1/2  in  continuum  flow 
(ref.  16).  Levin  (ref.  8)  has  shown  that  F  should  assume  values  of  0.  5  to  1. 1  for  condi¬ 
tions  of  meteor  flight  most  likely  to  be  encountered.  Conditions  for  typical  photographic 
meteors  are  best  represented  by  r=  1. 1  (free  molecular  flow)  (ref.  10).  For  the  larger 
asteroidal  meteors,  the  conditions  of  continuum  flow  prevail  with  r=  0.  5. 

The  dimensionless  shape  factor  A,  as  defined  in  equation  (2),  is  a  function  of  the 
shape  of  the  meteoroid.  For  a  sphere,  A  =  1.  208,  and  for  a  cube  traveling  with  one  face 
normal  to  the  direction  of  motion,  A  =  1.  00.  For  an  oblate  spheroid  with  the  minor  axis 
in  the  direction  of  flight,  ♦ 


where  b  is  the  ratio  of  the  minor  to  the  major  axis.  Realistic  estimates  of  A  based  on 
experimental  data  are  impossible  to  obtain.  Probable  limiting  values  of  A  are  difficult 
to  determine  since  a  meteoroid  can  be  conceived  to  have  a  shape  ranging  from  a  long 
rodlike  structure  to  a  platelike  structure. 

Experimentally,  Cook,  Jacchia,  and  McCrosky  (ref.  15)  found  that  if  FA  =  0.  92,  the 
average  luminosity  coefficient  of  the  three  asteroidal  meteors  would  be  in  agreement  with 
the  Trailblazer  (artificial  meteor)  experiment.  The  values  of  the  luminous  efficiency  for 
the  three  asteroidal  particles  plotted  in  figure  2  were  calculated  by  assuming  each  was  a 
sphere  (A  =  1.  208).  An  average  value  for  the  luminosity  coefficient  for  the  three  meteors 
was  calculated  based  on  the  linear  dependence  of  t  on  v  in  equation  (19),  and  this  value 
was  then  compared  with  the  luminosity  coefficient  derived  from  the  Trailblazer  experi¬ 
ment.  By  adjusting  the  drag- shape  factor  of  the  average  asteroidal  meteor  to  0.  92,  exact 
agreement  between  the  derived  luminosity  coefficients  of  Trailblazer  and  the  asteroidal 
meteor  average  was  found,  and  the  average  value  for  A  was  1.  84.  This  value  implies 
that  the  ’’average”  shape  of  the  associated  asteroidal  meteors  treated  was  an  oblate 
spheroid  with  the  major  axis  about  1.  9  times  the  minor  axis. 

In  treating  a  large  sample  of  cometary  meteor  data,  Verniani  (ref.  10)  indicated  that 
the  assumption  of  the  meteoroid  as  a  sphere  (A  =  1.  208)  is  a  good  estimate  since  the 
shape-factor  effect  will  average  out  if  the  meteoroid  shapes  are  assumed  to  be  entirely 
random.  ^  Individual  calculations  of  meteoroid  density  and  luminous  efficiency  are, 
therefore,  only  as  good  as  the  estimated  value  of  the  shape  factor.  For  a  large  sample 

^This  assumption  may  represent  a  lower  limit  for  A,  since  a  sphere  has  the  least 
ratio  of  mass  to  surface  area  of  any  solid  object. 
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treated  as  a  whole,  however,  the  effects  of  values  of  A  differing  from  that  of  a  sphere 
are  assumed  to  "average”  out,  but  for  individual  cases,  gross  errors  in  the  calculated 
results  are  possible. 


METEOROID  CHARACTERISTICS 
Density 

Until  recently,  the  average  density  of  meteoroids  was  assumed  to  be  the  same  as  the 
average  density  of  meteorites.  This  assumption  led  the  astronomer  to  believe  meteor¬ 
oids  were  stony  iron  or  iron  in  nature  with  a  density  in  the  range  of  3.  5  to  8.  0  grams  per 
cubic  centimeter.  With  this  assumption  of  meteoroid  densities,  investigators  in  the  field 
of  meteor  astronomy  were  able  to  calculate  atmospheric  densities  from  photographic  me¬ 
teor  data.  With  the  onset  of  rocketry  came  direct  measurement  of  the  atmospheric  den¬ 
sity  profile  at  altitudes  previously  calculated  from  meteor  data.  Comparison  of  the  cal¬ 
culated  and  measured  values  of  atmospheric  density  showed  that  the  calculated  values 
were  consistently  low  in  the  altitude  range  above  7  5  kilometers.  This  discrepancy  of  val¬ 
ues  could  be  resolved  if  the  density  of  meteoroids  were  assumed  to  be  much  lower  than 
that  previously  used.  Whipple  (ref.  18),  through  his  icy -conglomerate  comet  model,  pro¬ 
posed  just  such  a  low-density  fragile  structure  for  an  average  meteoroid  of  cometary  ori¬ 
gin.  Opik  was  also  a  proponent  of  the  fragile- structure  school  of  thought.  His  concept 
differs  slightly  from  that  of  Whipple  in  that  he  proposed  that  the  cometary  particles  enter¬ 
ing  the  Earth's  atmosphere  and  producing  meteors  are  dust  balls  composed  of  a  great 
number  of  individual  grains  (refs.  19  and  20).  In  effect,  both  models  imply  a  low  overall 
bulk  density  for  a  cometary  meteoroid. 

Whipple  (ref.  13)  and  others  then  reversed  the  process  of  calculation  by  assuming  at¬ 
mospheric  densities  to  be  known  and  meteor  densities  to  be  unknown.  This  procedure, 
with  improved  estimates  of  the  constants  involved  in  the  physical  equations  of  meteor  phe¬ 
nomena  and  recognition  of  meteor  fragmentation,  is  the  one  used  today.  Reevaluation  of 
meteoroid  densities  based  on  this  reversal  in  calculation  led  to  the  proposed  low  bulk 
density  of  a  meteoroid. 

Whipple,  in  his  latest  work  (ref.  17),  has  done  this  for  a  group  of  meteors  previously 
reduced  by  Jacchia  (ref.  12).  As  explained  previously,  the  right  side  of  equation  (8)  is  a 
function  only  of  the  photographic  data  and  atmospheric  density.  Whipple  refers  to  this 

_3 

value  as  K  .  Thus,  from  equation  (8), 

^A  meteoroid  is  a  particle  in  space;  a  meteorite  is  a  meteoroid  that  has  survived 
the  passage  through  the  Earth's  atmosphere  and  has  impacted  the  Earth.  Meteorites  are 
generally  believed  to  be  asteroidal  in  origin. 
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(21) 


From  the  data  of  reference  12,  Whipple  obtained  an  average  value  of  the-right  side  of  equa- 

fi  nr? 

tion  (21)  based  on  88  meteor  cases.  The  value  of  K  =  10  •  (units  of  zero  magnitude 
(visual)  g~^cm~^sec^)  together  with  the  values  of  =  10  and  FA  =  0.  92  led 

Whipple  to  conclude  that  the  typical  photographic  meteor  density  is  0.  44  gram  per  cubic 
centimeter. 

There  seems  to  be  an  inconsistency  in  the  values  and  therefore  in  the  proposed  den¬ 
sity  of  reference  17.  The  value  of  ==  10"  is  a  logarithmic  average  of  the  value 

of  Tq  determined  from  the  asteroidal  particles  in  reference  15  and  the  value  of  de¬ 
termined  from  the  Trailblazer  experiment  (ref.  14).  These  values  were  determined  with 
the  assumption  of  FA  =  0.  6  (corresponding  to  a  sphere  with  F  =  1/2).  As  explained  pre¬ 
viously,  if  FA  =  0.  92  is  assumed  for  the  asteroidal  particles,  the  values  of  for  the 
Trailblazer  and  the  asteroidal  particles  could  be  made  to  agree.  It  therefore  seems  that 
if  the  logarithmic  average  of  is  employed,  the  value  of  FA  =  0.  6  must  necessarily 
be  used,  or  if  FA  =  0.  92  is  to  be  used,  the  value  of  for  the  asteroidal  particles  and 
the  value  of  for  the  Trailblazer  experiment  must  agree  with  one  another  and  no  aver¬ 
aging  is  necessary.  Apparently  Whipple  recognized  this  inconsistency  because,  on  page  29 
of  reference  17  he  states  that  a  better  value  of  is  10"  •  ,  with  FA  =  0.92.  These 

latter  values  are  consistent  with  the  results  of  reference  15.  Employing  these  latter  val¬ 
ues,  Whipple  revised  his  best  estimate  of  meteoroid  density  to  0.  25  gram  per  cubic  centi¬ 
meter. 

Verniani  (refs.  10  and  21)  recently  made  an  exhaustive  study  of  a  total  of  700  meteor 
cases  in  order  to  define  meteoroid  densities.  In  reference  10,  equation  (8)  is  applied  to 
unpublished  data  obtained  by  Jacchia.  The  359  cases  treated  in  this  reference  had  mea¬ 
sured  decelerations  that  made  the  use  of  equation  (8)  possible.  The  meteors  were  cor¬ 
rected  for  fragmentation,  and  the  atmospheric  density  for  application  in  the  right  side  of 
equation  (8)  was  determined  from  the  known  height  of  the  meteor.  Since  each  meteor  was 
assumed  to  be  a  sphere,  FA  =  1. 3 ,  and  Verniani’s  derived  value  of  =  10“  ’  was 

used. 

The  results  of  Verniani's  calculations  are  shown  in  table  I.  These  average  values 
are  believed  to  be  significant  because  of  the  application  of  the  spherical  shape  factor, 
which,  as  explained  previously,  is  assumed  to  be  a  valid  assumption  when  applied  to  a 
large  sample.  Table  I  groups  the  meteors  into  subsamples  of  all  sporadic  meteors, 
sporadic  meteors  with  the  aphelion  distance  Q  less  than  7  astronomical  units,  sporadic 
meteors  with  the  aphelion  distance  greater  than  7  astronomical  units,  and  the  various 
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TABLE  I.  -  DENSITY  OF  METEORS  COMPUTED  BY  USE  OF 


DRAG  EQUATION  (REF.  21)  FROM  UNPUBLISHED  DATA 


Meteor 

Number 

Density, 
g/cu  cm 

All  sporadic  meteors 

247 

0.21 

Sporadic  meteors  with  aphelion  distance  Q  <  7  AU 

155 

.23 

Sporadic  meteors  with  aphelion  distance  Q  >  7  AU 

92 

.  16 

Taurids 

23 

.20 

Geminids 

20 

.92 

Aquarids 

16 

.24 

a  Capricornids 

13 

.  12 

Quadrantids 

9 

.  15 

Perseids 

9 

.28 

Orionids 

8 

.  14 

T]  Cygnids 

4 

.  14 

Lyrids 

3 

.28 

a  Hydrids 

3 

.30 

Virginids 

2 

-.70 

Draconids 

2 

-.001 

Total 

359 

identifiable  showers.  The  separation  of  the  sporadic  meteors  into  these  two  groups  is  an 
attempt  to  show  the  differences  between  long-  and  short-period  meteors  by  defining  arbi¬ 
trarily  an  aphelion  distance  of  7  astronomical  units  as  the  dividing  line.  The  results 
show  that  the  density  is  generally  smaller  than  0.  44  gram  per  cubic  centimeter,  with  the 
average  of  all  the  sporadic  meteors  slightly  greater  than  0.  2  gram  per  cubic  centimeter. 
This  value  agrees  with  the  results  as  finally  derived  in  reference  17  of  an  average  den¬ 
sity  of  0.  25  gram  per  cubic  centimeter. 

The  short-period  meteors  have  a  density  about  20  percent  greater  than  the  long- 
period  meteors.  Also,  the  density  of  shower  meteors  with  small  orbits  (Geminids, 
Taurids,  Quadrantids,  and  Aquarids)  is,  on  the  average,  larger  than  that  of  showers 
with  elongated  orbits  (Orionids  and  Perseids).  The  significance  of  this  difference  be¬ 
tween  long-  and  short-period  meteor  densities  is  uncertain.  Some  investigators  believe 
that  the  short-period  meteors  have  as  their  origin  the  inner  core  of  the  comet.  The 
short-period  comet  is,  on  the  average,  closer  to  the  Sun  than  the  long-period  comet  and 
will  disintegrate  more  rapidly.  The  long-period  meteors,  however,  originate  from  the 
outer  surfaces  of  the  long-period  comet  and  may  be  less  compact. 

A  search  among  the  meteors  studied  in  this  sample  revealed  no  dependence  of  the 
density  on  mass  or  photographic  magnitude  (ref.  21).  Furthermore,  the  meteors  used 
by  Whipple  (ref.  17)  from  which  he  obtained  a  density  of  0.  25  gram  per  cubic  centimeter 
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TABLE  11.  -  DENSITY  OF  METEORS  COMPUTED  BY  APPROXIMATE 
GROUP  METHOD  OF  CALCULATION  (REF.  21)  FROM 
PHOTOGRAPHIC  DATA  OF  REFERENCE  24 


Meteor 

Number 

Density, 
g/cu  cm 

All  sporadic  meteors 

280 

0.22 

Sporadic  meteors  with  aphelion  distance  Q  <  7  AU 

203 

.23 

Sporadic  meteors  with  aphelion  distance  Q>  7  AU 

77 

.  19 

Taurids 

17 

.  17 

Geminids 

16 

.32 

Aquarids 

10 

.49 

Orionids 

12 

.  18 

Perseids 

9 

.  14 

Quadrant  ids 

5 

.27 

Bielids 

2 

-.35 

a  Capricornids 

2 

-.  10 

Monocerotids 

1 

-.20 

Draconids 

1 

-.001 

Total 

355* 

are  of  a  much  greater  average  mass  than  those  used  in  reference  21,  which  produced  a 

density  of  0.  2  gram  per  cubic  centimeter.  This  range  of  mass  in  which  essentially  no 

“2 

dependence  of  density  on  mass  was  found  is  approximately  from  100  to  10  grams,  or 
four  orders  of  magnitude. 

As  pointed  out  previously,  it  is  not  always  possible  to  determine  meteoroid  decelera¬ 
tion  from  photographic  plates.  This  difficulty  often  makes  using  equation  (8)  impossible 
for  determining  meteoroid  density.  In  reference  21  an  approximate  method  is  developed 
for  determining  meteoroid  densities  that  does  not  require  meteor  deceleration  but  rather 
meteor  height  at  maximum  light  intensity.  The  method  relates  the  meteor's  shape, 
mass,  and  density  to  the  atmospheric  pressure  at  maximum  light  by  making  certain  as¬ 
sumptions  as  to  the  physical  strength  of  the  meteor.  It  is  indicated  that  the  method  is 
approximate  and  valid  only  when  applied  to  a  large  group  of  meteors.  The  method,  as 
derived,  gives  average  density  directly  and  loses  meaning  when  applied  to  single  cases. 

For  comparison  purposes,  table  II  lists  values  of  average  density  for  a  group  of  355 
meteors  (from  ref.  21)  calculated  by  this  method.  The  values  calculated  by  this  method 
are  in  close  agreement  with  the  direct  calculations  shown  in  table  I  for  the  sporadic  me¬ 
teors.  As  is  to  be  expected,  the  showers  calculated  by  the  two  different  methods  do  not 
show  good  agreement  with  one  another  when  the  sample  size  is  small.  These  results  in¬ 
dicate  an  average  value  of  meteoroid  density  typical  for  sporadic  meteors  to  be  about 
0.  2  gram  per  cubic  centimeter. 
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structure 


As  mentioned  in  the  preceding  section,  meteors,  irrespective  of  size,  were  earlier 
thought  to  be  compact  stones  or  iron  chunks.  With  better  knowledge  of  the  upper  atmo¬ 
sphere  and  observations  of  visual  and  photographic  meteors,  it  was  recognized  that  all 
meteoroids  could  not  have  such  a  high  density.  Reinterpretation  of  the  photographic  me¬ 
teor  data,  as  explained  in  the  preceding  section,  has  led  to  the  belief  that  meteoroids 
may  have  a  bulk  density  of  about  0.  2  gram  per  cubic  centimeter.  Furthermore,  photo¬ 
graphic  data  indicate  that  most  meteors  break  up  at  the  high  altitudes  at  which  they  first 
appear.  Opik  (ref.  11)  indicates  that  the  visible  meteors  typically  have  a  crushing 
strength  of  about  10^  dynes  per  square  centimeter  (0.  145  Ib/sq  in.),  or  about  the  crush¬ 
ing  strength  of  cigar  ash.  The  structure  and  composition  of  meteoroids  are  therefore 
also  important  considerations. 

The  spectra  of  many  meteors  have  been  obtained  (ref.  4),  but  no  specific  percentage 
compositions  are  reported.  The  occurrence  of  the  lines  of  specific  elements  in  the  me¬ 
teor  spectra  analyzed,  however,  is  given.  The  presence  of  the  following  elements  has 
been  detected:  sodium,  potassium,  silicon,  aluminum,  magnesium,  calcium,  manga¬ 
nese,  chromium,  and  iron  (ref.  22).  Calcium  and  iron  have  been  recognized  as  the  two 
most  prevalent  elements  in  meteor  spectra.  These  data,  however,  do  not  allow  any  con¬ 
clusions  regarding  the  typical  cometary  composition  because  of  their  qualitative  nature. 
In  order  to  obtain  quantitative  meteor  spectra,  the  luminous  efficiency  of  each  element 
in  the  spectrum  must  be  known. 

McCrosky  and  Soberman  (ref.  14)  have  indicated  another  means  for  estimating  the 
composition  of  typical  cometary  meteoroids.  The  relative  ratio  of  heavy  elements  to 
oxygen  in  stony  meteorites  is  approximately  the  same  as  the  relative  ratio  of  heavy  ele¬ 
ments  to  oxygen  in  the  Sun.  Since  the  meteor  spectra  indicate  the  presence  of  these 
heavier  elements,  the  assumption  is  made  that  the  same  relative  ratio  of  heavy  elements 
to  oxygen  exists  in  the  typical  cometary  meteoroid  as  in  the  Sun  and  in  stony  meteorites 
(ref.  14).  For  this  reason  it  has  been  concluded  that  the  typical  meteor  is  "stony”  in 
composition. 

The  observed  data  indicate  then  that  typically  a  meteor  appears  at  high  altitudes 
(from  which  low  density  and  low  crushing  strength  is  inferred),  and  the  assumption  is 
made  that  its  composition  is  typically  stony.  The  meteoroid  structure  proposed  must 
necessarily  produce  phenomena  in  agreement  with  the  observed  data  and  assumptions. 

Two  possible  explanations  exist  for  the  observed  phenomena:  the  first  is  drag  dis¬ 
tortion  and  the  second  is  the  "dust-ball"  concept  of  a  meteoroid.  Under  the  drag- 
distortion  interpretation,  a  compact  stony  meteoroid  entering  the  atmosphere  melts  be¬ 
cause  of  friction  into  a  spherical  drop.  Atmospheric  drag  causes  the  meteor  to  flatten 
out  into  an  oblate  spheroid  with  its  minor  axis  in  the  direction  of  flight.  The  increased 


ablation  cross  section  (frontal  area)  would  then  cause  the  meteor  to  appear  at  extremely 
low  atmospheric  densities.  Opik  (ref.  11),  however,  shows  that  for  compact  meteors 
under  the  influence  of  drag  distortion,  the  length  of  the  meteor  trail  should  decrease  with 
increasing  brightness.  No  such  effect  has  been  observed  statistically,  and  he  therefore 
concludes  that  the  explanation  of  the  increased  area  of  ablation  as  an  effect  of  deforma¬ 
tion  of  a  drop  is  not  acceptable.  If,  then,  the  compact  meteor  model  deforming  under  drag 
distortion  is  discarded,  the  acceptance  of  low  bulk  density,  low  crushing  strength,  and 
stony  composition  as  typical  of  photographic  meteors  is  necessary.  It  is  easy  to  visual¬ 
ize  this  type  of  particle  as  being  constructed  like  a  cigar  ash  with  layers  of  soft,  loosely 
bound  material,  or  as  a  stony  sponge. 

These  cigar  ash  or  stony  spongelike  particles  could  have  the  comets  as  their  origin 
as  proposed  in  the  icy- conglomerate  comet  model  of  reference  18.  In  essence,  the  pro¬ 
posal  is  that  the  comet  can  be  thought  of  as  a  large  frozen-gas  snowball  interlaced  with 
either  the  stony  sponge  or  cigar  ash  material  presented  previously.  Upon  repetitive  pas¬ 
sage  in  proximity  to  the  Sun,  the  gases  boil  off  and  leave  the  porous,  fragile,  low  bulk 
density  material  on  the  outer  surface  of  the  comet.  Solar  radiation,  internal  pressures, 
and  thermal  shocks  cause  the  outer  surface  of  the  comet  to  break  up  into  the  meteoroids 
observed  as  meteors.  This  particular  model  of  meteoroid  formation  seems  to  propose 
a  particle  that  satisfies  the  observed  phenomena. 

As  mentioned  previously,  Opik  also  proposed  a  meteoroid  model  to  meet  the  re¬ 
quirements  of  low  density  and  low  crushing  strength  (refs.  19  and  20).  He  proposed  a 
model  in  which  the  meteoroid  is  composed  of  a  great  number  of  individual  compact  grains 
loosely  held  together;  however,  he  offers  no  detailed  explanation  of  the  injection  mecha¬ 
nisms  of  such  a  particle.  The  differences  between  Opik's  dust-ball  meteoroids  and 
Whipple’s  cometary  particles  are  insignificant,  and  one  could  justify  the  same  type  of 
origin  for  the  dust  ball.  In  effect,  the  three  meteoroid  structures  previously  discussed 
could  be  grouped  as  one:  cometary  dust  balls. 

An  interesting  proposition  has  been  offered  by  Bonn  (ref.  23)  for  the  origin  of  mete¬ 
oroids  that  would  give  them  a  somewhat  different  type  of  construction.  In  that  reference, 
it  is  surmised  that  original  meteoroid  grains  were  condensed  from  the  primordial  gas  of 
the  solar  nebula  when  the  planets  were  formed.  The  grains  then  agglomerated  with  a  vol¬ 
atile  binding  material  to  form  comets.  The  method  of  injection  of  these  cometary  grains 
into  meteor itic  particles  then  is  the  same  as  that  proposed  by  the  icy- conglomerate 
comet  model. 

The  significant  contributions  of  the  theory  of  reference  23  is  in  the  proposed  struc¬ 
ture  of  the  cometary  meteoroid.  It  is  shown  that  the  particles  formed  by  primordial  gas 
condensation  would  grow  as  whiskers,  filaments,  platelets,  or  some  combination  result¬ 
ing  in  a  complex,  irregular  structure.  The  original  particle  can  then  be  described  as  a 
steel-wool  type  of  structure  as  one  possibility  or  a  many-pointed  starlike  structure  as 
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another.  The  platelet  structure  would  again  propose  the  ashlike  conglomerate.  Very 
little  change  in  the  original  particle  structure  would  result  from  cometary  agglomeration, 
and  thus  the  cometary  particles  of  today  would  have  very  nearly  preserved  the  primordial 
condensate  structure  and  composition  they  originally  possessed  (ref.  23).  This  theory  of 
meteoroid  structure  would  again  fulfill  the  observed  low  crushing  strength  and  low  den¬ 
sity  typical  of  photographic  meteoroid  particles. 

All  the  previously  discussed  possible  structures  for  typical  cometary  meteoroids 
have  two  characteristics  in  common:  distributed  mass  with  low  bulk  density  and  low 
crushing  strength.  The  question  then  arises  as  to  the  impact  damage  capabilities  of  this 
type  of  particle.  Existing  experimental  hypervelocity  impact  data  indicate  a  very  weak 

relation  between  penetration  and  particle  density  (p  cc  g_  ^  These  exist¬ 

ing  data,  however,  were  obtained  with  solid  homogeneous  projectiles  and  the  density- 
damage  relation  for  heterogeneous,  low  bulk  density  type  particles  is  unknown.  Further¬ 
more,  the  effect  of  the  low  crushing  strength  of  the  meteoroid  on  penetration  damage  is 
also  unknown.  Again,  all  experimental  and  theoretical  penetration  data  to  date  have  been 
obtained  with  projectiles,  the  crushing  strength  of  which  is  orders  of  magnitude  higher 
than  the  speculated  crushing  strength  of  the  typical  cometary  meteoroid.  In  general,  low 
bulk  density  and  low  crushing  strength  in  a  heterogeneous  particle  are  codependent,  and 
thus  the  individual  effects  are  not  likely  to  be  separated. 


Velocity 

A  knowledge  of  meteoroid  velocities,  besides  being  required  for  meteor  data  reduc¬ 
tion  purposes,  is  necessary  as  an  input  for  design  protection  calculations  (ref.  1). 
Ideally,  a  spatial  velocity  distribution  may  exist  for  meteoroids.  To  date,  a  spatial  ve¬ 
locity  distribution  curve  has  not  been  developed.  For  this  reason,  ’’average"  meteoroid 
velocities  have  been  proposed.  It  must  be  remembered,  however,  that  this  average  ve¬ 
locity  proposal  is  the  velocity  relative  to  the  Earth  and  not  necessarily  the  "impact"  ve¬ 
locity,  or  velocity  relative  to  the  space  vehicle  in  question.  For  these  reasons  a  discus¬ 
sion  of  meteoroid  average  velocities  is  felt  necessary. 

Figure  3(a)  presents  a  velocity  histogram  of  285  sporadic  meteor  cases  as  reported 
by  Hawkins  and  Southworth  (ref.  24).  The  velocities  are  relative  to  the  Earth  at  the 
upper  limit  of  the  atmosphere.  These  data  are  significant  to  study  since,  as  discussed 
in  a  following  section,  they  are  the  data  used  by  Whipple  in  proposing  his  latest  influx 
rates  (ref.  17) .  The  distribution  of  these  sporadic  data  is  bimodal  and  shows  a  peak  at  a 
velocity  of  about  20  kilometers  per  second  and  another  at  about  67  kilometers  per  second. 
The  reduced  values  between  45  and  55  kilometers  per  second  represent  meteors  in  rela¬ 
tively  unlikely  orbits  as  explained  in  reference  24. 
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A  velocity  histogram  derived  from  the  data  of  McCrosky  and  Posen  (ref.  25)  is  shown 
in  figure  3(b).  The  data  consist  of  2048  meteor  cases  that  were  determined  to  be  sporad¬ 
ic.  The  virtual  agreement  of  the  two  histograms  thus  verifies  the  randomness  of  the  data. 
As  indicated  in  reference  25,  the  method  of  reduction  used  for  these  latter  data  is  not  so 
accurate  as  the  method  used  in  reference  24.  The  randomness  of  the  errors,  however, 
validates  the  treatment  of  the  sample  as  a  group. 

An  arithmetic  average  of  the  two  histograms  yields  approximately  30  kilometers  per 
second  for  both  the  data  of  reference  24  and  those  of  reference  25.  This  value  of  average 
meteoroid  velocity  (relative  to  the  Earth)  has  frequently  been  used  for  the  impact  velocity 
in  meteoroid  protection  calculations. 

The  average  velocity  obtained  from  figure  3  represents  an  average  velocity  of  the  ob¬ 
served  sample.  The  observed  velocity  distribution,  however,  is  not  considered  to  be 
representative  of  the  true  velocity  distribution  in  space  because  of  a  velocity  bias  of  ob¬ 
servation.  Such  a  bias  exists  because,  in  general,  it  is  easier  to  ’’see”  the  high-velocity 
meteors  because  they  produce  brighter  trails;  that  is,  relatively  more  of  the  existing  me¬ 
teors  are  observed  when  their  velocities  are  high. 

As  indicated  in  reference  25,  the  velocity  bias  of  observation  on  a  photographic  plate 
is  the  combination  of  two  effects.  The  first  effect  is  seen  from  equation  (4),  with 
T  =  T^v.  This  relation  indicates  that  the  light  intensity  per  unit  time  from  a  meteor  of  a 
particular  mass  is  a  function  of  the  cube  of  its  velocity.  Comparison  of  light  energies 
emitted  by  individual  meteors  proportional  to  the  ratio  of  the  cube  of  their  velocities  is 
exact,  however,  only  when  the  masses  and  duration  times  are  equal  (i.  e. ,  when  dm/dt 
is  constant). 

The  second  factor  contributing  to  the  velocity  bias  of  observation  arises  from  the  ef¬ 
fective  exposure  time  on  the  photographic  plate.  The  plate  exposure  time  is  an  inverse 
function  of  the  translational  velocity  of  the  point  source;  therefore,  the  intensity  of  light 
recorded  on  the  plate  is  also  inversely  proportional  to  the  source  velocity.  This 
exposure-time  effect  will  influence  the  observed  velocity  histograms  in  a  manner  oppo¬ 
site  to  that  of  the  previous  velocity- energy  effect,  since  the  faster  the  meteor  travels, 
the  less  it  is  exposed  on  the  photographic  plate.  The  velocity  histograms  are  thus 
weighted  in  the  direction  of  the  low-velocity  meteors.  In  reality,  however,  the  situation 
is  not  quite  so  simple,  since  the  meteor  does  not  emit  light  energy  as  a  traveling  point 
source.  In  theory,  the  visible  light  energy  is  emitted  by  excited  meteoric  mass  parti¬ 
cles  caused  by  high-energy  collisions.  These  excited  particles  trail  the  meteoroid  by 
some  finite  distance  and  emit  as  a  large  number  of  point  sources.  The  exact  dependence 
of  exposure  time  on  meteor  velocity  is  further  complicated  by  the  existence  of  a  direc¬ 
tional  variation  of  the  velocity  vector  of  the  meteor  with  respect  to  the  photographic  plate 
and  by  the  variation  in  meteor  range.  (High-velocity  meteors  have,  statistically,  a 
greater  range. ) 
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The  exact  observational  bias  function  can  be  calculated  only  for  individual  cases  and 
is  quite  involved.  An  approximate  velocity  bias  for  large  samples,  however,  may  be 
chosen  as  proportional  to  1/v  by  combining  the  preceding  two  effects  as  recommended 
in  reference  25.  The  observed  count  is  then  adjusted  for  this  effect.  In  accepting  this  re¬ 
lation  for  the  removal  of  the  velocity  bias  of  the  histograms  of  figure  3,  it  must  be  as- 
sumed  that,  in  the  photographic  mass  range  of  interest  (10  to  10  g),  there  is  no  pro¬ 
nounced  dependence  of  meteoroid  velocity  on  mass.  The  further  assumption  is  made  that 
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these  histograms  are  representative  of  any  small  mass  range  within  the  limits  of  10  to 
10*^  gram.  Then,  for  this  "constant-mass”  histogram  there  is  a  higher  probability  of 
observing  the  high-velocity  particles  than  the  low-velocity  particles. 

In  establishing  a  more  nearly  correct  count,  it  is  assumed  that  all  existing  particles 
are  observed  for  velocities  between  70  and  75  kilometers  per  second  in  figure  3,  and  a 
correction  for  the  velocity  bias  of  the  photographic  observations  proportional  to  1/v  is 
made  for  the  other  velocity  ranges.  On  this  basis,  the  "true”  distribution  of  particle 
velocities  for  the  samples  considered  assume  the  shapes  shown  in  figure  4.  As  expected, 
these  histograms  indicate  many  "unobserved”  particles  of  low  velocity,  and  significantly 
this  indicates  that  there  is  more  flux  present  than  observed.  The  average  velocity  of  the 
true  distribution  of  the  data  of  reference  24  (fig.  4(a))  is  20.  2  kilometers  per  second  and 
of  the  data  of  reference  25  (fig.  4(b))  is  20.  3  kilometers  per  second,  values  that  are  sig¬ 
nificantly  less  than  the  nominal  30  kilometers  per  second  of  the  ’’observed”  distributions. 


MASS  INFLUX  RATE 

The  design  engineer  must  know  how  many  meteoroids  of  a  particular  mass  will  be 
encountered  by  the  vehicle  during  its  trajectory.  The  determination  of  meteoroid  influx 
rates  from  photographic  meteor  data  has  been  thoroughly  treated  by  Hawkins  and  Upton 
(ref.  26).  The  classical  approach  employed  in  reference  26  is  generally  accepted  and 
will  be  briefly  explained  herein.  The  photographic  data  obtained,  as  explained  previ¬ 
ously,  yield  the  light  intensity  and  photographic  magnitude  of  the  meteor  directly.  An 
examination  of  these  data  shows  an  increase  in  the  number  of  the  faint  meteors  as  com¬ 
pared  with  the  bright  meteors,  which  indicates  an  increase  in  the  number  of  meteors 
with  decreasing  mass. 

Hawkins  and  Upton  (ref.  26),  from  a  random  sample  of  285  sporadic  meteor  trails, 
derived  an  equation  for  the  cumulative  influx  rate  N,  expressed  as  a  function  of  the 
photographic  magnitude  over  the  range  of  photographic  magnitudes  from  -2.  0  to 
+4.  0,  in  the  form 


log  N  =  aM  -  log  Nq 

ir  ^ 


(22) 
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where  the  constant  a  is  the  slope  of  the  curve  of  log  N  against  Mp  and  Nq  is  the  num¬ 
ber  of  meteors  with  magnitudes  equal  to  or  less  than  zero.  Here,  the  photographic  mag¬ 
nitude  M  is  an  average  value  over  the  meteor  trail.  By  simply  counting  the  number  of 
meteors  with  magnitudes  less  than  or  equal  to  zero,  log  Nq  was  determined: 

log  Nq  =  -4.  34  (23) 

p  1 

(where  Nq  is  in  km"  hr"  ).  This  number  was  based  on  an  average  rate  of  observation 
by  the  viewing  cameras  of  0.  278  meteor  per  hour  over  a  collection  area  of  5980  square 
kilometers  of  sky.  The  number  of  meteors  increased  by  a  factor  of  3.  4  per  magnitude 
and  thus  produced  a  value  of  a  of  0.  537.  Equation  (22)  then  becomes 

log  N  =  0.  537  M  -  4.  34  (24) 

Jr 

The  representation  of  meteor  flux  rate  in  terms  of  photographic  magnitude  given  by 
equation  (24)  is  perhaps  the  most  accurate  form  of  this  type  of  information.  It  is  uncon¬ 
taminated  by  further  calculations  and  is  strictly  a  derived  representation  of  raw  data. 

The  accuracy  of  the  equation  is  a  function  only  of  the  accuracy  of  measurement  of  the 
photographic  magnitude  and  the  area  of  the  sky  photographed.  This  representation  of  me¬ 
teoroid  influx  rate,  however,  is  not  very  usable  for  space-vehicle-protection  require¬ 
ments,  since  a  knowledge  of  meteoroid  mass-influx  rates  is  necessary. 


Hawkins  and  Upton  Rate 

Hawkins  and  Upton  (ref.  26)  indicate  that  they  also  calculated  individual  masses  for 
each  sporadic  meteor  in  their  sample  of  285  sporadic  meteors.  This  calculation  can  be 
made,  as  explained  previously,  by  the  use  of  equation  (7)  with  an  assumed  value  of  the 
luminosity  coefficient.  The  choice  of  a  particular  value  for  the  luminosity  coefficient 
then  implies  that  a  consistent  value  of  an  average  meteoroid  density  is  determined.  This 
value  is  approximately  3  grams  per  cubic  centimeter  for  the  influx  rate  of  Hawkins  and 
Upton.  They  found  that  the  cumulative  influx  rate  N  could  be  expressed  as  a  function  of 
the  meteoroid  mass  m  in  the  form 

log  N  =  a^^  log  m  -  log  Nj^  (25) 

where  the  constant  a^  is  the  slope  of  the  curve  of  log  N  against  log  m  and  Nj^  is  the 
number  of  meteors  with  masses  equal  to  or  greater  than  1  gram.  Determination  of  the 
constants  a^  and  log  Nj^  yielded 
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log  N  =  -1.  34  log  m  -  2.64 


(26) 


-  2  - 1 

(where  N  is  in  km  hr  ).  Equation  (26)  is  a  more  useful  form  of  the  meteoroid  influx 
representation  for  a  space-vehicle  designer.  It  is  not  so  accurate  a  representation  of  the 
data  as  equation  (24),  since  it  does  have  implicit  in  it  an  assumed  luminosity  coefficient. 


Whipple  1963A  Rate 


Equations  (25)  and  (26)  are  the  basic  equations  used  by  Whipple  (ref.  17)  in  determin¬ 
ing  his  proposed  mass-influx  rate,  referred  to  as  the  "Whipple  1963A  flux  rate.  "  The 
basic  difference  between  this  proposed  influx  and  equation  (26)  as  used  by  Hawkins  and 
Upton  is  the  luminosity  coefficient  used  in  calculating  meteoroid  masses.  Since  this  co¬ 
efficient  is  the  basic  difference,  Whipple  adopted  equation  (26)  and  corrected  it  to  a  new 
luminous  efficiency  in  the  following  manner:  Equation  (7),  with  t  =  t^v,  shows  the  mass 
m  as 


(7) 


Since  the  masses  of  Hawkins  and  Upton  (ref.  26)  differ  from  the  mass  of  Whipple  (ref.  17) 
only  by  a  different  luminosity  coefficient,  the  ratio  between  the  two  mass  values  is 


”^W  _^o,H 
“h  ^o,W 


or 


o,W/ 


(27) 


where  the  subscripts 


Substituting  log  mjj 


W  and  H  refer  to  Whipple  and  Hawkins,  respectively. 

^o  H 

log  m^  =  log  — 2 —  +  log  mjj 

'^o,W 

from  equation  (28)  into  equation  (26)  yields 


Then 

(28) 
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log  N  =  1.  34  log 


-  1.  34  log  m  -  2.  64 


(29) 


T 


T 


o,H 


o,W 


which  can  be  used  to  correct  the  Hawkins  and  Upton  equation  to  any  value  of  the  luminos¬ 
ity  coefficient. 

^  PO  Q 1  “19  42 

Hawkins  and  Upton  used  =  10"  .  Whipple  (ref.  17)  proposed  =  10  '  ; 

equation  (29)  therefore  becomes,  for  the  Whipple  values, 

log  N  = -1.  34  log  m  -  4.  63  (30a) 

(where  N  is  in  km  hr  )  or 

log  N  = -1.  34  log  m  -  14.  18  (30b) 

(where  N  is  in  m"^sec“  V  Equations  (30a)  and  (30b)  (differing  only  by  their  units),  give 
the  proposed  influx  rate  of  Whipple  (ref.  17),  the  Whipple  1963A  flux  rate.  The  relations 
of  equations  (30)  do  not  contain  the  Earth  shielding  factor  of  1/2  that  Whipple  included  in 
his  proposed  influx  rate.  The  reduction  in  flux  by  a  factor  of  1/2  as  proposed  by  Whipple 
(ref.  17)  is  not  realistic  for  the  general  case,  since  obviously  the  effect  of  Earth  shield¬ 
ing  is  a  function  of  the  altitude  of  the  vehicle,  and  can  be  1/2  only  when  the  vehicle  is  at 
the  effective  meteoroid  boundary  of  the  Earth's  atmosphere. 

Included  in  the  proposed  influx  rate  of  equations  (30)  is  an  average  meteoroid  den¬ 
sity  consistent  with  the  assumed  luminosity  coefficient.  This  value  is  0.  44  gram  per 
cubic  centimeter.  The  influx  rate  of  equations  (30),  with  an  average  velocity  of  30  kilo¬ 
meters  per  second  and  an  average  density  of  0.  44  gram  per  cubic  centimeter,  has  be¬ 
come  the  most  widely  used  meteor  flux  model  to  date. 


Revised  Rate 

As  discussed  in  three  previous  sections,  better  values  of  meteoroid  mean  density 
(based  on  a  new  luminosity  coefficient)  and  meteoroid  average  velocity  appear  to  be  avail¬ 
able.  It  seems  reasonable,  therefore,  to  propose  a  meteoroid  flux  model  based  on  these 
values.  In  general,  equation  (26)  can  be  corrected  to  a  new  value  of  the  luminosity  coef¬ 
ficient  Tq  by  using 

,f^-20.91 

log  N  =  - 1.  34  log  m  -  2.  64  +  1.  34  log -  (31) 
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Verniani  (ref.  10)  arrived  at  =  10  ‘  and  =  0.  2  gram  per  cubic  centi¬ 

meter  as  representative  of  cometary  meteors.  If  this  value  of  is  accepted  as  the 
best  estimate  to  date,  equation  (31)  is  corrected  to  yield 

log  N  =  -1.  34  log  m  -  4.  23  (32) 

(where  N  is  in  km  hr  ).  As  shown  previously,  the  average  velocity  of  the  observed 
sample  (fig.  3,  p.  21)  is  about  30  kilometers  per  second.  The  use  of  equation  (32)  with 
an  average  velocity  of  30  kilometers  per  second  and  an  average  density  of  0.  2  gram  per 
cubic  centimeter  would  yield  a  consistent  flux  model  for  the  observed  sample. 

As  discussed  previously,  the  velocity  bias  of  observation  indicates  that  many  parti¬ 
cles  are  unobserved.  Since  the  design  engineer  is  more  interested  in  the  actual  flux  dis¬ 
tribution  than  in  the  observed  distribution,  equation  (32)  should  be  further  corrected  to 
account  for  these  inferred  "unobserved”  particles.  The  mass  influx  rate  derived  by 
Hawkins  and  Upton  (ref.  26)  has  in  it  a  correction  of  observation  based  on  the  assumption 
that  all  the  large  mass  particles  are  observed.  The  following  observational  bias  removal 
assumes  that  all  the  high-velocity  meteors  were  observed,  and  that  regardless  of  mass 
the  low-velocity  meteors  went  largely  unobserved. 

The  velocity- corrected  distribution  yielded  v^  =  20  kilometers  per  second  compared 
with  the  observed  distribution  value  of  of  30  kilometers  per  second.  The  decrease 

in  average  velocity  was  attributed  to  the  inferred  increase  in  the  number  of  low-velocity 
particles.  This  increase  in  particle  number  was  assumed  previously  to  be  proportional 
to  1/v  for  constant  particle  mass.  The  increase  in  sample  number,  for  essentially 
constant  mass,  would  thus  be  proportional  to  the  square  of  the  ratio  of  the  sample  aver¬ 
age  velocities;  that  is. 


or 


A  log  N  =  0.  352 


(33) 


Equation  (32)  is  then  modified  to  yield 

log  N= -1.34  log  m  -  3.88  (34) 

-2  - 1 

(where  N  is  in  km  hr  ).  Equation  (34)  with  =  0.  2  gram  per  cubic  centimeter  and 
V  =  20  kilometers  per  second  is  then  a  consistent  meteoroid  flux  model  for  a  near-Earth 
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Cumulative  influx  rate,  N,  km'^hr' 


Mass,  m,  g 


Figure  5.  -  Deduced  influx  rates  based  on  data  of  Hawkins  and  Upton  (ref.  26). 


environment  representing  an  estimate  of  the  '•actual*'  rather  than  the  observed  distribu- 

-2  0 

tion.  The  mass  range  over  which  this  model  can  be  assumed  to  apply  is  10  to  10  gram, 
which  is  the  range  of  the  observed  data. 

The  various  flux  models  discussed  herein  along  with  their  associated  values  of  den¬ 
sity  and  velocity  are  shown  in  figure  5.  The  flux  models  of  Hawkins  and  Upton,  Whipple 
196 3A,  and  the  solid  curve  of  equation  (32)  are  all  representative  of  the  observed  refer¬ 
ence  data,  and  the  differences  are  wholly  attributed  to  the  assumed  luminosity  coefficient. 
The  dashed  curve  is  the  inferred  "actual"  distribution  as  discussed  previously  in  equa¬ 
tion  (34). 

The  revised  actual  distribution  (dashed  curve)  indicates  an  increase  in  the  meteor 
hazard  as  compared  with  the  Whipple  1963A  curve.  The  proposed  associated  lower  den¬ 
sity  and  velocity,  however,  counteract  much  of  this  increase,  as  shown  in  a  later  section. 


EARTH  SHIELDING  FACTOR 

When  a  flux  model  is  used  for  meteoroid  protection  calculations  for  an  Earth  satel¬ 
lite,  the  influence  of  the  Earth  as  a  shield  should  be  considered.  Whipple  recognized 
this  influence  but  applied  a  single  factor  of  1/2.  In  reality,  the  Earth  shielding  factor  E 
is  a  variable  and  can  be  so  computed  as  a  function  of  the  altitude  of  the  satellite  above 
the  Earth’s  surface  that,  at  any  altitude, 

=  EN  (35) 

o 

where  N  is  obtained  from  equation  (31). 

In  calculating  E,  it  will  be  assumed  that  the  meteoroid  population  is  isotropic.  Fig¬ 
ure  6  represents  a  unit  hemisphere  of  space  of  unit  radius,  with  its  center  at  the  satel¬ 
lite;  r  is  the  true  radius  of  the  Earth,  h  is  the  satellite  altitude  above  the  surface  of  the 
Earth,  and  R  is  the  "effective"  radius  of  the  Earth.  Using  the  effective  radius  of  the 
Earth  recognizes  that  most  meteors  burn  out  in  the  upper  atmosphere  at  altitudes  of  70 
to  100  kilometers.  The  atmosphere  of  the  Earth  therefore  provides  an  effective  shield 
to  meteoroids  and  thereby  serves  to  increase  the  radius  of  the  occluding  body.  For 
practical  purposes,  the  atmosphere  at  90  kilometers,  or  approximately  50  nautical  miles, 
should  provide  an  effective  shield  to  meteoroid  passage;  that  is,  Ar  =  50  miles.  With  the 
use  of  this  assumption,  the  effective  Earth  shielding  as  a  function  of  altitude  can  be  de¬ 
rived  from  reference  to  the  geometry  of  figure  6. 

The  shielding  factor  of  a  satellite  at  altitude  h  above  the  Earth's  surface  is  deter¬ 
mined  to  be 
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(36) 


E  -  1  ,  V(r  +  h)^  -  (r  + 

2  2(r  +  h) 

Figure  7  shows  the  Earth  shielding  factor  as  a  function  of  altitude  based  on  the  as¬ 
sumption  that  the  atmosphere  is  an  effective  shield  up  to  50  nautical  miles.  The  Earth's 
occlusion  is  of  importance  to  satellites  in  orbit  at  less  than  about  2500  nautical  miles. 


DESIGN  IMPLICATIONS 


From  the  preceding  discussion  it  is  obvious  that  the  reduction  of  photographic  meteor 
data  into  meaningful  information  useful  to  the  space-vehicle  designer  is  not  a  simple  and 
precise  task.  Furthermore,  there  are  many  limitations  and  assumptions  involved  in  the 
reduction  methods  that  can  have  notable  effects  and  implications  on  the  design  protection 
requirements.  These  limitations  and  assumptions  can  best  be  appreciated  by  examining 
a  typical  relation  used  for  predicting  required  armor  thickness  to  prevent  critical  mete¬ 
oroid  damage. 

Reference  1  presents  an  equation  for  required  single- sheet  armor  thickness  as  a 
function  of  mission  parameters  and  meteoroid  hazard  as 


\PtJ  VC/ 


-In  P(0) 


1/3/3 


1/3/3 


+  2y 


In  terms  of  only  meteoroid  properties, 
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(38) 


where  and 
factors  a  and 


V  are  the  meteor  bulk  density  and  impact  velocity,  respectively. 
-)3  are  related  to  the  meteoroid  flux  according  to  the  equation 


The 


N  =  Q!m~^ 


(39) 


The  exponent  -/3  is  the  slope  of  the  general  curve  of  log  N  against  log  m  (eq.  (25))  and 
a  is  the  same  as  Nj^. 
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Meteoroid  Density 


In  a  previous  section  the  calculation  of  the  density  for  a  particular  meteoroid  was 
shown  to  be  dependent  on  several  factors  other  than  the  data  obtained  from  the  photo¬ 
graphic  plate.  In  particular,  the  luminosity  coefficient  and  the  drag-area  coefficient 
rA,  as  well  as  the  correct  application  of  the  fragmentation  index,  were  shown  to  be  highly 
influencing  factors.  The  uncertainty  introduced  in  the  value  of  the  meteoroid  density  by 
the  uncertainty  in  these  factors  is  easily  seen  from  equation  (11)  as  proportional  to 

VtqCTA^.  It  is  impossible,  however,  to  place  uncertainty  limits  on  either  or  TA, 
since  at  best  they  are  educated  estimates.  The  treatment  of  large  groups  of  data,  as  ex¬ 
plained  previously,  is  assumed  to  nullify  the  uncertainty  in  the  drag-area  coefficient  and 
thus  to  reduce  the  uncertainty  in  the  "average"  meteoroid  density  primarily  to  a  function 
of  Tq  (if  fragmentation  is  assumed  to  be  correctly  applied).  Fortunately,  equation  (38) 
shows  that  the  required  armor  thickness  for  space-vehicle  protection  bears  a  very  weak 
relation  to  the  meteoroid  density.  This  weak  relation,  which  was  derived  for  solid  homo¬ 
geneous  particles,  may  be  a  conservative  relation  for  the  low-bulk- density,  low-crushing- 
strength  meteoroid  particle. 


Meteoroid  Velocity 

The  treatment  of  estimated  impact  velocity  warrants  further  study.  The  attempt 
herein  to  remove  the  observational  bias  of  the  data  has  resulted  in  a  reduction  of  the 
average  velocity  of  the  meteor  sample  relative  to  the  Earth  from  a  nominal  30  to  20  kilo¬ 
meters  per  second.  The  correction  presented  is  not  precise  but  rather  is  an  attempt  to 
show  an  effect  that  can  markedly  influence  the  total  meteoroid  hazard.  If  an  average  me¬ 
teoroid  velocity  relative  to  the  Earth  is  to  be  better  defined,  a  more  precise  treatment  is 
needed  in  which  the  observational  bias  is  removed  from  individual  meteor  cases  by  con¬ 
sidering  all  the  factors  involved  and  the  fact  that  all  meteors  are  not  of  the  same  mass. 
The  treatment  presented  assumed  that  the  velocity  histograms  were  independent  of  mass. 
It  is  probable  that  large  particles  will  be  observed  regardless  of  their  velocity;  this  pos¬ 
sibility  makes  the  correction  presented  invalid  for  this  region  of  mass. 

The  value  of  20  kilometers  per  second  deduced  herein  as  a  typical  or  average  mete¬ 
oroid  velocity  is  relative  to  the  Earth  at  the  limit  of  its  upper  atmosphere  and  is  not  nec¬ 
essarily  the  value  that  should  be  used  in  equation  (38).  Ideally,  the  value  in  equation  (38) 
should  be  the  effective  meteoroid  impact  velocity  relative  to  the  space  vehicle  for  its 
mission.  Uncertainties  in  relative  impact  velocity  affect  the  required  single- sheet  armor 
thickness  drastically  (proportional  to  v^^^),  and  a  significant  effort  is  therefore  indi¬ 
cated  for  the  determination  of  accurate  values.  Meteoroid  relative  impact  velocity  is  a 
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function  of  the  vehicle  trajectory  in  space  and  the  spatial  meteoroid  velocity  distribution. 
The  second  of  these  two  items  is  at  present  undefined.  Much  work  remains  in  the  defini¬ 
tion  of  this  spatial  velocity  distribution.  The  heliocentric  as  well  as  the  geocentric  rela¬ 
tions  involved  are  necessary,  and  a  method  of  incorporation  of  the  results  into  protection 
equations  such  as  equation  (38)  must  be  devised  for  a  precise  analysis  of  the  hazard.  To 
date,  meteor  data  are  lacking  for  the  definition  of  this  distribution.  For  an  Earth- orbiting 
vehicle  at  high  altitude  (long  periods),  a  20-kilometer-per- second  relative  impact  velocity 
for  particles  in  the  photographic  mass  range  may  be  an  adequate  estimate.  More  advanced 
missions,  however,  such  as  those  that  leave  an  Earth  orbit  or  that  involve  orbits  with 
short  periods  or  high  eccentricity,  will  require  a  more  sophisticated  analysis  of  the  esti¬ 
mation  of  meteoroid  impact  velocities. 

Mass-Influx  Rate 

The  influence  of  the  mass-influx  rate  on  the  required  armor  thickness  in  equation  (38) 
is  indicated  by  the  values  of  a  and  /3. 

The  value  of  jS,  the  slope  of  the  curve  of  log  N  against  log  m,  equal  to  1.  34  estab¬ 
lished  by  Hawkins  and  Upton  in  1958  (ref.  26)  has  remained  unchanged.  Although  a  more 
precise  value  of  /3  can  probably  be  determined  from  an  increased  sample  mass  range 
and  sample  number,  jS  =  1.  34  can  be  assumed  to  be  sufficiently  valid  for  use  over  the 
photographic  mass  range. 

Uncertainties  in  a  will  reflect  uncertainties  in  6  approximately  proportional  to  the 
1/4  power  (if  jS  =  1.  34  is  assumed).  The  value  of  a  is  seen  from  equation  (31)  to  be  di¬ 
rectly  related  to  the  uncertainty  in  the  luminosity  coefficient  raised  to  the  1.  34  power. 
Comparing  the  influx  rates  of  Whipple  (eq.  (30))  and  Hawkins  and  Upton  (eq.  (26))  shows 
a  two-order-of-magnitude  difference  in  the  value  of  a  (from  2.  64  to  4.  63)  and  corre¬ 
spondingly  in  the  value  of  flux.  This  variation  is  wholly  attributable  to  the  differences  in 
the  values  used  for  t^.  This  indicated  relation  reflects  the  strong  dependence  of  the  re¬ 
quired  armor  thickness  on  the  absolute  value  of  t^.  As  explained  previously,  this  abso¬ 
lute  value  is  difficult  to  determine;  therefore,  confidence  in  (and  finally  a)  becomes 
related  to  the  "educated  guess”  made  in  determining  t^. 

In  addition  to  uncertainties  in  a  (related  to  the  accuracy  of  t^),  the  value  of  the  flux 
will  depend  on  the  magnitude  of  the  velocity  bias  of  the  observation  and  the  effective  re¬ 
moval  of  this  bias  from  the  observed  sample.  The  dashed  curve  of  figure  5  (p.  29)  indi¬ 
cates  an  "actual"  distribution  inferred  from  observational  data  based  on  an  approximate 
removal  of  the  observational  bias,  as  explained  previously.  This  representation  of  an 
"actual"  rather  than  the  observed  influx  rate  is  currently  felt  to  be  more  meaningful  for 
protection  calculations.  A  better  curve  could  be  inferred  from  a  larger  sample  of  data 


by  removal  of  the  bias  from  individual  meteor  cases.  The  ’’probability”  of  observation, 
therefore,  does  not  have  to  be  made  independent  of  mass,  as  explained  previously. 

Finally,  the  influx  rates  discussed  herein  should  specifically  be  limited  in  use  to  a 
near- Earth  environment.  The  effect  of  gravity  focusing  by  the  Earth  is  implicit  in  the 
measured  flux  rate  and  therefore  in  the  inferred  actual  flux  rate.  Since  the  focusing  ef¬ 
fect  tends  to  concentrate  meteoroid  particles  in  the  immediate  vicinity  of  the  Earth,  flxix 
rates  could  be  significantly  different  from  the  observed  values  at  large  distances  from  the 
Earth.  The  focusing  effect,  being  velocity  sensitive,  requires  the  treatment  of  individual 
meteor  cases  for  greater  precision  in  the  determination  of  the  flux. 


Comparison  with  Whipple  1963A  Hazard 


The  meteoroid  hazard  as  proposed  herein  based  on  the  inferred  ’’actual”  flux 
and  revised  values  of  density  and  velocity  can  be  compared  with  the  Whipple  1963A  hazard 
on  the  basis  of  equation  (38).  The  ratio  of  the  required  armor  thickness  for  the  revised 
hazard  based  on  the  dashed- curve  environment  model  proposed  herein  to  the  armor 
thickness  given  by  the  Whipple  1963A  hazard  model  6^  for  the  same  mission  is  given  by 
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(40) 


The  following  table  lists  the  various  input  values  for  each  model: 


Model 

Meteoroid 

density, 

Pm’ 

g/cu  cm 

Average 

meteoroid 

velocity, 

V, 

km/sec 

Slope  of  mass 
against  influx 
rate, 

P 

log  a 

Whipple  1963A 

0.44 

30 

1.34 

-4.63 

Revised  environment 

.2 

20 

1.34 

-3.88 

(dashed  curve,  fig.  5) 

Substituting  into  equation  (40)  the  various  values  of  p^,  v,  /3,  and  a  for  each  model 
yields 


Sr 

-^=1.03  (41) 

5w 
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The  revised  hazard  indicates  that  about  3  percent  more  armor  is  needed  than  predicted  by 
the  Whipple  1963A  hazard  model.  This  minor  increase  in  required  armor  thickness  is 
due  to  the  counterbalancing  effects  of  an  increased  indicated  flux  and  a  decreased  average 
impact  velocity.  For  all  practical  purposes  the  two  models  yield  the  same  armor  re¬ 
quirements. 


Meteoroid  Structure  and  Shape 

Two  major  factors  have  arisen  in  the  determination  of  the  average  meteoroid  density 
that  can  have  a  further  marked  influence  on  the  required  protection.  These  factors  re¬ 
late  to  the  shape  and  the  structure  of  the  meteoroids  and  the  basis  of  available  impact  re¬ 
lations.  Current  knowledge  of  penetration  and  damage  effects  due  to  the  impact  of  hyper¬ 
velocity  particles  (e.  g. ,  eq.  (38))  has  been  obtained  from  experimental  and  theoretical 
studies  based  on  solid  homogeneous  particles,  usually  spherical  or  cylindrical  (equal 
length  and  diameter)  in  shape,  and  with  relatively  large  crushing  strengths.  Current  es¬ 
timates  of  the  meteoroid  properties,  as  described  herein,  however,  indicate  that  mete¬ 
oroid  particles  may  be  substantially  different  from  the  man-made  particles  in  both  struc¬ 
ture  and  shape. 

First,  it  is  believed  that  meteoroids  may  vary  considerably  in  shape  and  range  possi¬ 
bly  from  long  needlelike  rods  to  flat  plates.  Experimentally,  the  shape  of  a  solid  impact¬ 
ing  particle  and  its  orientation  at  impact  have  been  shown  to  have  a  major  effect  on  the 
resulting  penetration  (e.g. ,  ref.  27).  If  these  same  effects  persist  for  meteoroids,  sig¬ 
nificant  variations  in  impact  damage  may  result  from  individual  meteoroid  particles  of 
the  same  nominal  characteristics,  that  is,  mass,  density,  and  velocity. 

Secondly,  and  perhaps  more  important,  meteoroid  particles,  by  virtue  of  their  known 
constituents,  low  bulk  density,  and  possible  origins,  are  believed  to  possess  relatively 
low  crushing  strength  with  a  structure  containing  a  heterogeneous  distribution  of  mass. 

If  this  is  indeed  true,  then  the  resulting  impact  effects  could  be  different  from  the  cur¬ 
rently  available  solid-particle  relations,  and,  hopefully,  less  severe.  Furthermore,  if 
such  differences  exist,  it  would  be  important  to  determine  the  possible  influence  on  these 
differences  of  such  factors  as  meteoroid  velocity,  shape,  and  mass.  Theoretical  and  ex¬ 
perimental  investigations  may  shed  some  light  on  possible  impact  differences,  while  ad¬ 
ditional  meteor  astronomy  techniques  or  spacecraft  experiments  may  be  useful  for  fur¬ 
ther  structure  and  shape  determinations. 


CONCLUDING  REMARKS 

It  is  shown  herein  how  photographic  meteor  data  and  methods  of  analysis  assembled 
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from  various  scientific  sources  are  evaluated  and  interpreted  for  use  in  the  design  and 
development  of  meteoroid  protection  for  large  spacecraft  components.  It  was  indicated, 
however,  that  considerable  uncertainties  remain  in  astronomically  derived  meteoroid 
properties.  Since  the  pertinent  meteoroid  properties  for  the  larger  particles,  which  con¬ 
stitute  a  significant  threat  to  space  devices  with  large  vulnerable  areas,  cannot  be  obtained 
by  practical  flight  experiments,  further  effort  should  be  directed  toward  the  reduction  of 
the  uncertainties  in  the  photographic  meteor  analysis  process.  Efforts  are  needed  to  in¬ 
vestigate  the  impact  characteristics  of  meteoroid- type  particles  as  discussed  herein  and 
to  define  more  accurately  meteoroid  structure,  shape,  and  composition. 

Lewis  Research  Center, 

National  Aeronautics  and  Space  Administration, 

Cleveland,  Ohio,  May  7,  1965. 
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APPENDIX  -  SYMBOLS 


A 

shape  factor 

Nq 

cumulative  influx  rate  of  zero- 

\ 

vehicle  vulnerable  area 

magnitude  meteors 

a 

slope  of  magnitude  influx  rate 

Nl 

cumulative  influx  rate  of  1-gram 

particles  (same  as  a) 

^1 

slope  of  mass  influx  rate 

1 

n 

exponent 

observed  meteor  acceleration 

ODS 

P 

penetration 

theoretical  meteor  acceleration 

t 

P(0) 

probability  of  no  critical  damage 

b 

ratio  of  minor  axis  of  meteoroid 

particle  to  major  axis 

Q 

aphelion  distance 

C 

sonic  velocity  of  impacted  target 

R 

effective  radius  of  Earth 

D 

constant  in  penetration  equation 

r 

mean  radius  of  Earth 

Dl 

constant 

Ar 

difference  between  effective  ra¬ 

1 

dius  and  mean  radius  of  Earth, 

E 

Earth  shielding  factor 

R  -  r 

h 

altitude  of  satellite  above 

S 

frontal  area  of  meteor 

Earth's  surface 

s 

mass  loss  parameter 

I 

intensity  of  light  given  off  by 

meteor 

®0 

mass  loss  parameter  at  particu¬ 

lar  point  on  meteor  trail 

r  /  /•  1 

-1/3 

K 

-(  1  /-idt 

T 

mission  time 

_  \pv2  dt;  J  ,3  _ 

t 

time 

“p 

photographic  magnitude 

te 

time  at  end  of  meteor 

m 

meteoroid  mass 

V 

meteoroid  velocity 

meteoroid  mass  at  end  of  meteor 

V 

average  meteoroid  velocity 

meteoroid  mass  before  it  enters 

"^c 

corrected  average  velocity 

Earth’s  atmosphere 

^obs 

observed  average  velocity 

N 

cumulative  influx  rate 

a 

cumulative  influx  rate  of  1-gram 

Nc 

corrected  cumulative  influx  rate 

particles  (same  as  Nj) 

Nobs 

observed  cumulative  influx  rate 

iS 

slope  of  mass  agsinst  influx  rate 

Ns 

Earth- shielded  cumulative  in- 

(1.34) 

flux  rate 

r 

particle  drag  coefficient 
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5 

e 


p 


required  armor  thickness  to  pre¬ 
vent  critical  damage 

constant  in  penetration  equation 
(2/3) 

atmospheric  density 
meteoroid  density 
target  material  density 


T  luminous  efficiency 

luminosity  coefficient 
X  fragmentation  index 

Subscripts: 


H 

Hawkins 

R 

revised 

W 

Whipple 

39 


REFERENCES 


1.  Loeffler,  I.  J, ;  Lieblein,  Seymour;  and  Clough,  Nestor:  Meteoroid  Protection  for 

Space  Radiators,  Progress  in  Astronautics  and  Aeronautics.  Vol.  E.  Academic 
Press,  Inc.,  1963,  pp.  551-579. 

2.  Dubin,  Maurice;  and  McCracken,  Curtis  W.  :  Measurements  of  Distribution  of  In¬ 

terplanetary  Dust.  Astron.  J. ,  vol.  67,  no.  5,  June  1962,  pp.  248-256. 

3.  Hastings,  Earl  C. ,  Jr. :  The  Explorer  XVI  Micrometeoroid  Satellite  -  Supplement  1, 

Preliminary  Results  for  the  Period  January  14,  1963,  through  March  2,  1963. 
NASATMX-824,  1963. 

4.  Watson,  F.  G.  :  Between  the  Planets.  Harvard  Univ.  Press,  1956. 

5.  Hawkins,  Geralds.:  The  Meteor  Population.  NASA  CR- 51365,  1963. 

6.  Whipple,  F.  L.  :  Particulate  Contents  of  Space.  Medical  and  Biological  Aspects  of 

the  Energies  of  Space,  P.  A.  Campbell,  ed. ,  Columbia  Univ.  Press,  1961. 

7.  Whipple,  Fred  L. ;  and  Jacchia,  Luigi  G.  :  Reduction  Methods  for  Photographic  Me¬ 

teor  Trails.  Smithsonian  Contributions  to  Astrophys. ,  vol.  1,  no.  2,  1957,  pp. 
183-206. 

8.  Levin,  Yuri.  V.  :  The  Physical  Theory  of  Meteors  and  Meteoroid  Matter  in  the  Solar 

System.  Pergamon  Press,  1963. 

9.  Jacchia,  L.  G.  :  The  Physical  Theory  of  Meteors.  VEI  -  Fragmentation  as  Cause  of 

the  Faint  Meteor  Anomaly.  Astrophys.  J. ,  vol.  121,  1955,  pp.  521-527. 

10.  Verniani,  Franco:  On  the  Luminous  Efficiency  of  Meteors.  NASA  CR- 55904,  1964. 

11.  Opik,  E.  J.  :  Physics  of  Meteor  Flight  in  the  Atmosphere.  Intersci.  Pub.,  Inc., 

1958. 

12.  Jacchia,  L.  G.  :  A  Comparative  Analysis  of  Atmospheric  Densities  from  Meteor  De¬ 

celerations  Observed  in  Massachusetts  and  New  Mexico.  Tech.  Rept.  No.  10, 
Harvard  College  Observatory,  1952. 

13.  Whipple,  F.  L.  :  On  Meteor  Masses  and  Densities.  Astron.  J. ,  vol.  57,  1956,  pp. 

28-29. 

14.  McCrosky,  R.  E. ;  and  Soberman,  R.  K.  :  Results  from  an  Artificial  Iron  Meteoroid 

at  10  km/sec.  Rept.  No.  AFCRL  X62-803,  Air  Force  Cambridge  Res.  Labs., 

July  1962. 

15.  Cook,  A.  F. ;  Jacchia,  L.  G. ;  and  McCrosky,  R.  C.  :  Luminous  Efficiency  of  Iron 

and  Stone  Asteroidal  Meteors.  Proc.  Symposium  on  Astronomy  and  Phys.  of  Me- 


40 


teors,  1963,  pp.  209-220. 

16.  McCrosky,  R.  E.  :  Observations  of  Simulated  Meteors.  Smithsonian  Contributions 

to  Astrophys. ,  vol.  5,  no.  4,  1961. 

17.  Whipple,  F.  L.  :  On  Meteoroids  and  Penetration.  Smithsonian  Astrophys.  Observa¬ 

tory,  1963. 

18.  Whipple,  F.  L.  :  A  Comet  Model.  I  -  The  Acceleration  of  Comet  Encke.  Astrophys. 

J. ,  vol.  Ill,  Mar.  1950,  pp.  375-394. 

19.  Opik,  E.  J.  :  The  Masses  and  Structures  of  Meteors.  Meteors,  T.  R.  Kaiser,  ed. , 

Pergamon  Press,  1955. 

20.  Opik,  E.  J.  :  Problems  in  the  Physics  of  Meteors.  Am.  J.  Phys. ,  vol.  26,  no.  2, 

Feb.  1958,  pp.  70-80. 

21.  Verniani,  F.  :  On  the  Density  of  Meteoroids.  II.  The  Density  of  Faint  Photographic 

Meteors.  Nuovo  Cimento,  ser.  10,  vol.  33,  no.  4,  Aug.  16,  1964,  pp.  1173-1184. 

22.  Dauvillier,  A.  :  Cosmic  Dust.  George  Newnes,  Ltd. ,  1963. 

23.  Donn,  Bertram:  The  Origin  and  Nature  of  Solid  Particles  in  Space.  Ann.  N.  Y. 

Acad.  Sci.,  vol.  119,  Nov.  11,  1964,  pp.  5-16. 

24.  Hawkins,  G.  S. ;  and  Southworth,  R.  B.  :  The  Statistics  of  Meteors  in  the  Earth’s 

Atmosphere.  Smithsonian  Contributions  to  Astrophys. ,  vol.  2,  no.  11,  1958,  pp. 
349-364. 

25.  McCrosky,  R.  E. ;  and  Posen,  A.  :  Orbital  Elements  of  Photographic  Meteors. 

Smithsonian  Contributions  to  Astrophys. ,  vol.  4,  no.  2,  1961,  pp.  15-84. 

26.  Hawkins,  G.  S. ;  and  Upton,  E.  K.  L.  ;  The  Influx  Rate  of  Meteors  in  Earth’s  Atmo¬ 

sphere.  Astrophys.  J. ,  vol.  128,  no.  3,  Nov.  1958,  pp.  727-735. 

27.  Charters,  A.  C.  :  A  Critique  of  Accelerator  Techniques  for  Hypervelocity  Impact. 

Proc.  Sixth  Symposium  on  Hypervelocity  Impact  (DDC  No.  AD  423798),  vol.  1, 

Aug.  1963,  pp.  1-40. 


NASA-Langley,  1965  E-2900 


41 


aeronautical  atid  space  activities  of  the  United  States  shall  he 
conducted  so  as  to  contribute  ...  to  the  expansion  of  human  knowl¬ 
edge  of  phenomena  in  the  atmosphere  and  space.  The  Administration 
shall  provide  for  the  widest  practicable  and  appropriate  dissemination 
of  information  concerning  its  activities  and  the  results  thereof.** 

— National  Aeronautics  and  Space  Act  of  1958 
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